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(54) ROAD SURFACE CONDITION ESTIMATION DEVICE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a road surface 
condition estimation device hard to be influenced by a 
noise. 

SOLUTION: A vehicle speed signal which is computed in 
a wheel speed calculating circuit 35 from a signal output 
from a wheel speed sensor 10 and which includes a 
vibration component at the time of impulse-excited 
wheel cylinder hydraulic pressure in an appointed cycle 
is input in a frequency computing circuit 30 and an 
amplitude computing circuit 32 for computing a 
frequency which is highly correlated with several wave 
form data, and an amplitude which corresponds to the 
frequency. An attenuation rate computing circuit 34 
computes an attenuation rate, and a road surface □ 
estimation device 14 judges a road surface □ gradient to 
be small when the attenuation rate is large and large 
when that is small on the basis of the attenuation rate 
computed in the attenuation rate computing circuit 34. 
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* NOTICES * 

JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1. This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2. **** shows the word which can not be translated. 
3.1n the drawings, any words are not translated. 

CLAIMS . 

[Claim(s)] , . ^ . 

[Claim 1] Correlation with a detection means and the data point of two or more predetermined 
frequencies from which the time series data and the frequency of a signal differ whenever 
[ wheel speed / in the condition that the oscillating component is contained in the signal 
whenever / wheel speed / with the 1st period ] is searched for whenever [ wheel speed / which 
detects whenever / wheel speed / and outputs a signal whenever / wheel speed ]. In the 
resonance frequency calculated with a resonance frequency operation means to calculate the 
resonance frequency of a data-point empty vehicle ring speed signal with high correlation, and 
said resonance frequency operation means It asks for the correlation coefficient of the time 
series data of a signal, and the data point of the frequency called for by said resonance 
frequency operation whenever [ wheel speed / in the condition that the oscillating component is 
contained in the signal whenever / wheel speed / the 2nd period shorter than said 1st period j. 
Road surface state estimation equipment including an amplitude operation means to calculate 
the amplitude of a signal whenever [ said wheel speed ] from the correlation coefficient for which 
it asked, and a road surface state estimation means to presume a road surface condition based 
on the degree of attenuation of the amplitude of a signal whenever [ wheel speed / which was 
calculated with said amplitude operation means ]. 

[Claim 2] Road surface state estimation equipment which includes further a wheel-cylinder oil 
pressure excitation means to excite wheel-cylinder oil pressure whenever [ said wheel speed J 
when the amplitude of a signal is small. 

[Translation done.] 
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JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1. This document has been translated by computer. So the translation may not reflect the original 
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2. **** shows the word which can not be translated. 
3.1n the drawings, any words are not translated. 



DETAILED DESCRIPTION 



[Detailed Description of the Invention] 

[0001] . . I . 

[Field of the Invention] This invention relates to road surface state estimation equipment, and .n 
case it presumes a road surface condition, for example, road surface mu inclinat.on, especially, it 
relates to robust road surface state estimation equipment to the oscillation frequency change by 
tire modification etc. by using two or more detection frequencies. 

[0002] . , 

[Description of the Prior Art] Vibration of whenever [ wheel speed ] changes to road surface mu 
inclination, and vibration of whenever [ wheel speed ] is continued comparat.vely long in the 
condition i.e., the condition that road surface mu inclination is large, that the wheel grips the 
road surface. On the other hand, in the small (near a grip limitation (i.e., road surface mu 
clination) of a wheel) condition, vibration of whenever [ wheel speed ] is converged 



in 



comparatively quickly. 

[0003] Therefore, if the degree of attenuation of signal empty vehicle ring rate vibration is 
measured whenever [ wheel speed ], a road surface condition, for example, road surface mu 
inclination, is detectable. , 
[0004] The degree of attenuation of vibration can be calculated [ whenever / wheel speed J from 
the amplitude value of a signal whenever [ wheel speed ]. However, if amplitude value is directly 
calculated from the instantaneous value of the data point of a signal whenever [ wheel speed J, 
since the sudden noise is contained in the data point, it is difficult to detect the degree of 
attenuation only with instantaneous value. 

[0005] It was made in order that this invention might cancel the above-mentioned trouble and it 
aims at offering ** and the road surface state estimation equipment which cannot be easily 
influenced of a noise using the technique of a correlation operation to detect the amplitude as a 
wave-like feature parameter. 

[0006] L . . 

[Means for Solving the Problem] In order to attain the above-mentioned purpose, this invention 
whenever [ wheel speed / which detects whenever / wheel speed / and outputs a signal 
whenever / wheel speed ] A detection means. Correlation with the data point of two or more 
predetermined frequencies from which the time series data and the frequency of a signal differ 
whenever [ wheel speed / in the condition that the oscillating component is contained in the 
signal whenever / wheel speed / with the 1st period ] is searched for. In the resonance 
frequency calculated with a resonance frequency operation means to calculate the resonance 
frequency of a data-point empty vehicle ring speed signal with high correlat.on, and said 
resonance frequency operation means It asks for the correlation coefficient of the time series 
data of a signal, and the data point of the frequency called for by said resonance frequency 
operation whenever [ wheel speed / in the condition that the oscillating component is contained 
in the signal whenever / wheel speed / the 2nd period shorter than said 1st period J. It 
constituted from a correlation coefficient for which it asked including an amplitude operation 
means to calculate the amplitude of a signal whenever [ said wheel speed ]. and a road surface 
state estimation means to presume a road surface condition based on the degree of attenuation 
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of the amplitude of a sigr^henever [ wheel speed / which was cWlated with said amplitude 
operation means ]. 

[0007] The vibration frequency operation means of this invention searches for correlation with 
two or more data points from which the time series data and vibration frequency of a signal 
differ whenever [ wheel speed / in the condition that the oscillating component is contained in 
the signal whenever / wheel speed / with the 1st period ], and calculates the resonance 
frequency of a data-point empty vehicle ring speed signal with high correlation. As this data 
point, the data of a sine wave and a cosine wave can be used. An amplitude operation means is 
the 2nd period shorter than the 1 st period, and calculates the amplitude of a signal in the 
vibration frequency calculated with the vibration frequency operation means whenever [ wheel 
speed / in the condition that the oscillating component is contained in the signal whenever / 
wheel speed ]. This amplitude can ask for the correlation coefficient of the time series data of a 
signal and the data point of the frequency called for by said frequency operation means 
whenever [ wheel speed ]. and can calculate it from the correlation coefficient for which it asked. 
As a data point of an amplitude operation means, the data of a different sine wave from the 
above-mentioned sine wave and a cosine wave and a cosine wave can be used. 
[0008] A road surface state estimation means presumes a road surface condition based on the 
degree, for example, the attenuation factor, or attenuation rate of attenuation of the amplitude 
calculated with the amplitude operation means. When the degree of attenuation is quick, it is 
easy to slide on a road surface condition, and since a wheel continues wheel vibration 
comparatively long in the condition of gripping the road surface and vibration is converged 
comparatively quickly near the grip limitation of a wheel as explained above, when the degree of 
attenuation is late, it can be judged that it is hard to slide on a road surface condition. 
[0009] It is effective, if a wheel-cylinder oil pressure excitation means to excite wheel-cylinder 
oil pressure is established further and wheel vibration is excited whenever [ wheel speed ] in this 
invention, when the amplitude of a signal is small. Thereby, an oscillating component comes to be 
contained in a signal whenever [ wheel speed ]. In addition, in this invention, wheel vibration may 
be produced using the disturbance from a road surface, without using a wheel-cylinder oil 
pressure excitation means. 
[0010] 

[Embodiment of the Invention] Hereafter, with reference to a drawing, the gestalt of operation of 
the road surface state estimation equipment of this invention is explained to a detail. 
[001 1] As shown in drawing 1 , the gestalt of this operation was attached so that a rotation 
gearing's gear tooth fixed to the wheel might be countered, and is equipped with the sensor 10 
whenever [ wheel speed / which outputs the electrical signal with which vibration frequency 
changes in proportion to whenever / wheel speed ]. Whenever [ this wheel speed ], although 
attached in each of four flowers of a car, the sensor 10 was represented and was illustrated one. 

[0012] Whenever [ wheel speed ], an arithmetic circuit 35 is connected to the outgoing end of a 
sensor 10 whenever [ wheel speed / which calculates whenever / wheel speed / from the period 
of vibration or the oscillating frequency of the vibrating electrical signal ], and an electrical signal 
or a numerical signal proportional to whenever [ wheel speed ] etc. is outputted to it The 
frequency arithmetic circuit 30 which calculates the frequency of a signal whenever [ wheel 
speed ] using the data point of two or more predetermined frequencies from which a frequency 
differs a long period, and the amplitude arithmetic circuit 32 which calculates the amplitude of a 
signal whenever [ to the frequency calculated the short period in the frequency arithmetic circuit 
30 / wheel speed ] are connected to the outgoing end of an arithmetic circuit 35 whenever 
[ wheel speed ]. This frequency arithmetic circuit 30 is connected so that the data of the 
calculated frequency may be inputted into the amplitude arithmetic circuit 32. 
[0013] The amplitude arithmetic circuit 32 is connected to the road surface mu inclination 
presumption equipment 14 which presumes road surface mu inclination from an attenuation 
factor through the attenuation factor arithmetic circuit 34 which calculates the attenuation 
factor of the amplitude. 

[0014] Furthermore, the output of the amplitude arithmetic circuit 32 is inputted into the brake 
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ic^Jltputs the control command to a hydrtBl 



hydraulic control 16 whicJ^Btputs the control command to a hydrlBE power package 18, when 
the amplitude in resonance frequency turns into amplitude below predetermined. Based on this 
control input, the hydraulic power package 18 equipped with the ABS control valve which 
switches the oil pressure of the master cylinder 26 according to the treading strength of a brake 
pedal 24 and the source of high pressure unrelated to treading strength which is not illustrated, 
and controls the oil pressure (wheel-cylinder oil pressure) of a wheel cylinder 20 is connected to 
the outgoing end of the brake hydraulic control 16. The ABS control valve consists of a boost 
bulb which boosts wheel-cylinder oil pressure, and a reducing valve which decompresses wheel- 
cylinder oil pressure. 

[0015] If the pressure sensor 22 which consisted of semiconductor pressure sensors for 
detecting wheel-cylinder oil pressure to a wheel cylinder 20 etc. can be attached, the property 
of attenuation factor-average wheel-cylinder oil pressure can also be used for presumption of a 
road surface condition. In this case, the output of a pressure sensor 22 is connected to road 
surface mu inclination presumption equipment 14. 

[0016] While being outputted from an arithmetic circuit 35 whenever [ wheel speed ] according 
to the gestalt of this operation, whenever [ containing the oscillating component when carrying 
out impulse excitation of the wheel-cylinder oil pressure a predetermined period / wheel speed ], 
a signal is inputted into the vibration frequency arithmetic circuit 30 and the amplitude arithmetic 
circuit 32, and the vibration frequency and the amplitude of a signal are the following whenever 
[ wheel speed ], and it is made and calculated. In addition, below, the sampling period of 
measurement is explained about 1 and the example which made the period of impulse excitation 
the 512 sample point (s), 0.512 [ i.e., ], (ms). 

[0017] First, an ABS control valve is controlled by the brake hydraulic control 16, and impulse 
excitation of the wheel-cylinder oil pressure is carried out a predetermined period with it. This 
impulse excitation is performed, only when it is the following and the amplitude of a signal is 
below a predetermined value whenever [ wheel speed / which was made and calculated ]. 
Drawing 2 shows change of whenever [ wheel-cylinder oil pressure / when it carries out impulse 
excitation the excitation period 0.512 (s), having used the average wheel-cylinder oil pressure 
command as 4 (MPa) /, and wheel speed ]. 

[0018] In addition, change of an oscillating component and an average oscillating component is 
shown in drawing 3 and drawing 4 whenever [ wheel speed / when carrying out impulse 
excitation of the average wheel-cylinder oil pressure command as 2, 4, 6, 8, 10, and 12 (MPa) ]. 
[0019] The operation of the resonance frequency of a signal is explained whenever [ by the 
vibration frequency arithmetic circuit 30 / wheel speed ]. Resonance frequency calculates 
correlation with the time series data of a signal whenever [ two or more sine wave and cosine 
wave, and wheel speed ], whenever [ wheel speed ], chooses the sine wave with highest time 
series data and correlation or cosine wave of a signal, and searches for it. [ from which the 
frequency defined beforehand differs ] 

[0020] First, correlation with the time series data of the excitation period 0.512 (s), and the sine 
wave of the wave number of 6, 7, 8, 9, and 10 and a cosine wave is calculated respectively. Sine 
wave C j m of the wave number m And cosine wave S j m It is expressed as follows. In addition, 
it may replace with the wave number and a sine wave and a cosine wave may be expressed using 
vibration frequency. 
[0021] 
[Equation 1] 

Cj m =(1/256) cos (2pim (2J+D/1024) ... (1) 
Sj m =(1/256) sin (2pim (2j+1)/1024) ... (2) 

j= 0, and 1 and 2 — if the time series data in a ... and 51 1i time are set to xn (n=i -51 1, i— 51 0, 
i) — real part of the correlation coefficient in i time Ri m And imaginary part Ii m It is given by 
the degree type. [ however, ] 
[0022] 
[Equation 2] 
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R i 1 — 2 Cd • B x. 

n ■ I -511 



I, 1 =2 So-' x, 

u-1 -511 



(3) 



(4) 



Sta l£ abrftha sum of aquaras Caquare of tha an^ituda) of a oorralation coaffician. 

is called for by the following formula. 

[0024] 

P i m =(R i m)2-Kl i m) 2 ... (5) 

The above (3) and (4) types can deform as follows. 

[0025] 
[Equation 3] 



i - i 

= £ C- 1 x. - Cr" x,- 6 ij + Cr" xi 

0-1-61J 

=R " +C ( Xl -Xl-512 ) 

1-1 

= E S.-' x„ - S,-" Xi-si, + S, ■■ X, 

». 1-612 

= 1 !-! » +S ( Xl -Xl-612 ) 



(6) 



(7) 



T00261 Since it can ask for the sum of squares of a correlation coefficient by the multiplication in 
Kd^erm, and addition with the 1st term and the 2nd term, it .can d™*. count of an 
operation sharply, so that I may be understood from the above (6) and 7) types 
[0027] The greatest sum of squares P i m among the sums of squares (m- 6. 7 8. 9. W P m ot 
the cor elatfon coefficient detected at the i time It chooses and is the greatest sum of squares 
Pirn The following formulas are followed from the wave number and it ,s the resonance 
frequency fm of a signal whenever [ wheel speed ]. It determines. 
[0028] 

ForTxlmpL 2 P Sm at the time of m= 8 If it is max, resonance frequency can be determined as fB 

[0 8 029°]' 5 A S 2 mlnl 2 ned H a Z bove. in order to calculate the correlation coefficient of the time series 
datf of a 572 sample point to the one wave number, resonance frequency wi.l ca.cu.ate a long 

lOOZO] In addition, resonance frequency can also be determined by the following formulas as a 
weighting average of P i m (m- 6, 7, 8, 9, 10). 
[0031] 
[Equation 4] 

10 io 

f d = (S P, 1 f . ) / (2 P. 1 ) 

Next in the resonance frequency of a signal, the amplitude of a signal is calculated [ whenever / 
wheel speed / which was detected the long period as mentioned above ] whenever wheel 
r P eed snort period in the amplitude arithmetic circuit 32. This amplitude ,s canted from a 
Correlation coefficient with the sine wave of vibration frequency and cos.ne wave wh.ch were 
called for from the above-mentioned vibration frequency arithmetic circu.t. 
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[0032] As for the vibratiij^equency of 6, 7, 8, 9, and 10, one per^Jlorresponds [ the wave 
number ] to the sample point of 85, 73, 64, 57, and 51 respectively to the time series data of the 
excitation period 0.512 (s). 

[0033] It is the sine wave c j m of the wave number m like detection of the vibration frequency 
in the long period explained above. And cosine wave s j m It expresses as follows. 
[0034] 
[Equation 5] 

cj m =(2-/Nm) cos (pi (2j+1)/Nm) ... (9) 
sj m =(2-/Nm) sin (pi (2j+1)/Nm) ... (10) 

j= 0. and 1 and 2 — it is .... Nm-1. Nm =85, and 73, 64, 57 and 51. [ however, ] 

[0035] Thus, the one wave number is received. Nm In order to calculate the correlation 

coefficient of a sample point, the amplitude will calculate a short period. 

[0036] And real part of the correlation coefficient in i time Ri m And imaginary part Ii m It 

calculates from a degree type. 

[0037] 

[Equation 6] 

i 

r t B = 2 c B - x. •••(11) 



* =2 s B ' x t • • • (1 2) 



however, n'=n mod Nm it is . 

[0038] From the above, the amplitude is called for by the following formula. 
[0039] 

ai m = {(r i m) 2+0 i m) 2} 2 ... (13) 
The above (1 1) and (12) types can deform as follows. 
[0040] 
[Equation 7] 
i-i 

r , Q = 2 c„ B x n - Ci- 1 xi-ju + d- m x, 

o*l -Hi 

=r f-i ' +c ,•■ ( x, -x,-„) ■ • • (14) 



i . 



I -1 

= 2 s n - 

□ - 1 -Si 



- Si- 1 XJ-8.+ Si * X| 



(1 5) 



however, i'=(i~N m) mod Nm =i mod Nm it is . 

[0041] By transforming a formula as mentioned above, as the case where vibration frequency 
was calculated a long period was explained, the count of an operation can be decreased sharply. 
[0042] ai m while it judges whether the amplitude is increasing in the attenuation factor 
arithmetic circuit 34 by comparing the amplitude in front of 1 period with this amplitude and the 
amplitude is increasing — > — ai-Nmm it is — since wheel-suspension vibration is continuing 
comparatively long, by the road surface input or the brake oil pressure control, to a case, it 
judges that vibration was added and the magnitude of attenuation ignores to it. 
[0043] It is an attenuation factor in order to, presume the road surface mu inclination as a road 
surface condition on the other hand, since it is being completed comparatively quickly by wheel- 
suspension vibration while the amplitude is decreasing, k i m is calculated as follows. 
[0044] the amplitude ai m, ai~1 m, and ... ai-Nmm The damping property to include can be 
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b^J^ng exponential functions like the cun^^c 



approximated with the foT^Jng exponential functions like the cunWiown in drawing 6 (A). 
[0045] 

y=exp sin (omega t+phi) (~ki m t) 

The envelope Y of this exponential function expresses amplitude fluctuation, and is expressed 
with the following formula. 
[0046] Y=exp(-ki m t) 

If the logarithm of the above-mentioned formula is taken, it is an attenuation factor ki m by 
searching for the inclination over t of logY, as it is set to logY=-ki m t and shown in drawing 6 
(B). It can ask. 

[0047] With road surface mu presumption equipment 14, a road surface condition is judged based 
on the attenuation factor calculated in the attenuation factor arithmetic circuit 34. When an 
attenuation factor is big (i.e., when vibration declines promptly), it is easy to slide on a road 
surface, i.e., road surface mu inclination is small, or since road surface mu inclination is small, a 
road surface mu judges that threshold value is approached. 

[0048] On the contrary, when an attenuation factor is small (i.e., when vibration continues for a 
long time), it is judged that it is hard to slide on a road surface, i.e., it is generous to the 
threshold value of a road surface mu since road surface mu inclination is road surface mu large 
or large. 

[0049] Furthermore, it can judge any of a low mu way and a quantity mu way transit road 
surfaces are by being able to calculate average wheel-cylinder oil pressure from the wheel- 
cylinder oil pressure measured, being able to create the table showing the damping property of 
suspension vibration of drawing 5 which shows the relation between average wheel-cylinder oil 
pressure and the damping factor calculated as mentioned above, and comparing this table, 
wheel-cylinder oil pressure at present, and a damping factor, if a pressure sensor 20 can be 
attached. 

[0050] Although the example which excites wheel-cylinder oil pressure by the ABS control valve 
was explained above, wheel-cylinder oil pressure is made to generate impulse-fluctuation, and 
w heel-suspension excitation may be made to cause using the actuator for TRC or VSC. 
[0051] Moreover, the amplitude may be extracted by removing average braking deceleration from 
a signal whenever [ wheel speed ], extracting only an oscillating component, starting a wave for 
every period of excitation, equalizing a wave in a rate region (for example, ten (rad/s) units), and 
searching for correlation with the sine wave and cosine wave which were explained above. 
[0052] The endurance of an actuator can be improved, while being able to fall the noise of an 
actuator according to the gestalt of this operation, without always exciting brake oil pressure 
since it is exciting only when the amplitude of vibration becomes small as explained above. 
Moreover, since the vibration frequency made into an object from a predetermined band is 
detected, road surface mu inclination can be presumed regardless of whether also when an 
oscillation characteristic changed with exchange of ignorance [ the oscillation characteristic of a 
wheel-suspension system ], or a tire etc., could respond, and TORAIBA operated the brake. 
[0053] In addition, except a smooth way, since wheel-suspension excitation can be made to 
cause by disturbance, the wheel-cylinder oil pressure excitation means explained with the 
gestalt of the above-mentioned implementation may be omitted. 

[0054] Moreover, although the gestalt of this operation explained wheel-suspension vibration 
(resonance frequency of about 20Hz) as an example, it is applicable not only to this but wheel- 
car-body vibration (resonance frequency of about 40Hz). 
[0055] 

[Effect of the Invention] Since two or more detection frequencies are used according to this 
invention in case a road surface condition, for example, mu inclination of a road surface, is 
presumed as explained above, the effectiveness that robust road surface state estimation 
equipment can be offered to the oscillation frequency change by tire modification etc. is 
acquired. 



[Translation done.] 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 



2005/01/24 



1/1 V 



* NOTICES * 




JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

LThis document has been translated by computer. So the translation may not reflect the original 
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DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[Drawing 1] It is the block diagram of the gestalt of operation of this invention. 
[Drawing 2] It is the diagram showing change of whenever [ wheel-cylinder oil pressure / when it 
carries out impulse excitation, having used the average wheel-cylinder oil pressure command as 
4 (MPa) /, and wheel speed ]. 

[Drawing 3] It is the diagram showing change of an oscillating component and an average 
oscillating component whenever [ wheel speed / when carrying out impulse excitation of the 
average wheel-cylinder oil pressure command as 2, 4, 6, and 8 (MPa) ]. 
[Drawing 4] It is the diagram showing change of an oscillating component and an average 
oscillating component whenever [ wheel speed / when carrying out impulse excitation of the 
average wheel-cylinder oil pressure command as 10 and 12 (MPa) ]. 

[Drawing 5] It is the diagram showing the table showing the damping property of a damping factor 
and suspension vibration. 

[Drawing 6] It is a diagram for explaining how to calculate the amplitude using an exponential 
function. 

[Description of Notations] 

10 It is Sensor whenever [ Wheel Speed ]. 

14 Road Surface Mu Inclination Presumption Equipment 

18 Hydraulic Power Package 

34 Attenuation Factor Arithmetic Circuit 



[Translation done.] 
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